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Abstract 

We present the result of the three loop anomalous dimension of non-singlet transversity 
operator in QCD for the Mellin moment N = 16. The obtained result coincides with the 
prediction from arXiv: 1203 . 1022 and can serve as a confirmation of the correctness of the 
general expression for three loop anomalous dimension of non-singlet transversity operator 
in QCD for the arbitrary Mellin moment. 



A general expression for the three loop anomalous dimension of non-singlet transverse 
operator in QCD for the arbitrary Mellin moment was obtained recently pQ. It has been 
reconstructed from the known first fifteen moments with the help of LLL-algorithm [2]. In 
spite of the fact that the reconstructed result seems to be correct it would be nice to have 
an additional confirmation from the direct calculations of higher moments, starting from 
N = 16. 

The local flavor non-singlet twist-2 transversity operator is given by 
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where = (z/2) [7^7" — 7 W 7 M ], A r are the Gell-Mann matrices for S"?7(3) flavoI ., cor- 
responds to the covariant derivative in QCD, q(q) denote the quark and antiquark fields, 
and the operator S symmetrizes the Lorentz indices and subtracts the trace terms. This 
operator appears, for example, in the study of semi-inclusive deeply inelastic scattering 
(SIDIS) 111 El E] and in the polarized Drell— Yan process El El M HQ • The scaling 
violation of the transversity distribution was explored in leading [HI [T21 EH EH E51 EE] and 
next-to- leading order [TTJ [181 EH]. At the three- loop order first 15 Mellin moments were 
obtained in Refs. [201 EU E21 E31 El E] • 

The calculation of the 16th moment for the anomalous dimension associated with 
transversity operator ([T]) is several times more difficult than the calculation of the corre- 
sponding moment for non-singlet anomalous dimension associated with the deeply inelastic 
structure functions F 2 (x, Q 2 ) and F L (x, Q 2 ), which was obtained at third order in Ref. [25] 
by method from Refs. [2E1 EZ] ■ To perform this kind of computation we had to improve 
substantially the method used for the calculation of lower moments in Ref. [TJ. In this 
paper we give a detailed description of our calculations. 

The anomalous dimension of the operator in MS-like scheme can be obtained from the 
corresponding renormalization constant by means of the following general formula 
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where g stands for strong gauge coupling, a is a gauge fixing parameter, and CrR ns (ot, g 2 ) 
is the coefficient of the lowest pole in e in the Laurent expansion of the renormalization 
constant for the operator 
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For calculation of the renormalization constants, following [28] (see also [291 EQl EI]), we 
use the multiplicative renormalizability of corresponding Green's functions. The renormal- 
ization constants relate the dimensionally regularized one-particle-irreducible Green's 
function with renormalized one as: 
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where g| and «b are the bare charge and the bare gauge fixing parameter, correspondingly, 
with 

oo 
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and Z 3 being the gluon field renormalization constant. 

In order to calculate the required anomalous dimension, following [32], we consider the 
Green's function G l ^J R ' ns , which is obtained by contracting the matrix element of the local 
operator ([!]) with the source term = A Ml . . . A™ 

u(p, s)G^ ns X r u(p, s) = J N (q t (p) | 0™^ „ x | gi(p)> , (6) 

where p and s denote the four-vectors of the momentum and spin of the external quark 
line, u(p, s) is the corresponding bi-spinor, A 2 = 0. The contraction with the source term 
J^v = A Ml . . . A MiV allows to write a general expression for the corresponding projector, 
which can be found in Ref. [52] . 

The unrenormalized Green's function has the following Lorentz structure 
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with unphysical constants c^| fc=1 5 . To determine £ TR ' ns (p) we use the following projector 
(see [23]) ' 

sTR ' nS W = - U Nc{A .p)^(D-2) { Tl[ ^ ^^ TR ' nS l " (A -p)Tr[p^-] 

+z(A-p)Tr[a w K^ TR ' nS ]}' ( 8 ) 

where iV c denotes the number of colors and D = 4 — 2e. 

Once the renormalization constant for E TR ' ns (p) is known, it is straightforward to com- 
pute Z TR,n S by means of 

TR,iis = Z ^ -Z^TR.ns, (9) 

with Z q being the renormalization constant for the quark field. 

As in our previous calculations [lj we made use of the program DIANA [33J, which calls 
QGRAF [3l] to generate all diagrams and the FORM package COLOR [35] for evaluation of 
the color traces. The unrenormalized three- loop £ TR ' ns (p) was computed with the help of 
the FORM [36] package MINCER [37J. Corresponding renormalization constant ZjjrR, u was 
determined from the requirement that the poles in e cancel in the r.h.s. of Eq. (j3J). 

Since the right-hand side of Eq. (j3J) contains the bare gauge fixing parameter qb, we 
should perform all calculations up to two loops with the arbitrary gauge fixing parame- 
ter a (i.e. the propagator of gluon is [g^ v — (1 — a)q ll q u /(l 2 )/ < l 2 ), while for the three-loop 
calculations it is sufficient to use the Feynman gauge a — 1. To obtain the result for the 
renormalization constant from Eq. we put a = 1 only after expansion in the right-hand 
side of Eq. g]). 
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It is obvious from the discussion presented above that our main task was to extract 
E TR ' ns (p) contribution from the unrenormalized Green function for the transersity operator. 
A convenient way to do this is to get rid of the contraction with A in Eqs. ([6]) and (EJ) 
and perform the Passarino — Veltman decomposition of the resulting tensor in the following 
way (see [32]): 
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where {} means the symmetrization with respect to the indices. It should be pointed 
out that Eq. (1101) has one additional index in comparison with the operator in Eq. (CQ) as 
projector (|SJ) itself contains an additional vector A. Due to the fact that the contraction 
with gives A 2 = it is sufficient to calculate only the coefficient A\ from Eq. ([TO]) . The 
corresponding projector reads [32] (note again, that the projector has 17 indices for the 
16th moment): 
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and the general expression for the arbitrary iV can be found in Ref. [32J. 

It turns out that for the calculation of the 16th moment the above-mentioned proce- 
dure has to be substantially optimized. The main difficulties for higher moments are the 
following. Firstly, the number of terms for an operator vertex with some outgoing gauge 
lines increases very fast (see Appendix of Ref. [32J). Secondly, the number of terms in the 
projectors grow even more rapidly. However, since the projector (|TT|) when applied to a 
diagram acts only on external Lorentz indices of the operator vertex, the contraction 

nOM.TK.Bs = n(J ,.., 17 0™« , fl 17 = V (14) 

is universal for all diagrams and we can substitute this expression (with proper relabeling) 
in place of the operator vertex. It should be stressed that such contraction produces only 
several Lorentz structures for arbitrary moment N. For example, in the case of an operator 
vertex without gauge lines attached, i. e., when all covariant derivatives in operator (JTJ) are 
substituted by the ordinary ones, we have 10 Lorentz structures multiplied by combinations 
of scalar products of internal and loop momenta 

no™-^(o) = Q ltx Q v Q) + Qi, P u n 2 (p, Q) + in^Q v n 3 (p, Q) + Pi^Pu tUp, Q) 

+ i^QQu Hs(p, Q) + IpQVv tUip, Q) + i^pQu Hiip, Q) + i^ppu H%{p, Q) 
+ i v nu'tUip,Q) +7*/)m^io(p><3) > (15) 
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Figure 1: Examples of diagrams containing three gauge fields in the operator vertex 



where Q corresponds to external momentum, and p stands for loop momenta flowing 
through the operator. In MINCER notations, p = pi, . . . ,p 8 . For operator vertex with one 
outgoing gauge line there are 72 Lorentz structures with 3 different indices, for operator 
vertex with two gauge fields we have 710 Lorentz structures with 4 different indices and for 
operator vertex with three gauge fields the number of structures with 5 different indices 
is 8900 (the total number is 9692). In general case the functions depend on scalar 
products of momenta entering to operator vertex. Thus, we have created a file containing 
the expressions for coefficients and substituted them in a given diagram only after 
contraction of Lorentz indices and computation of traces of 7-matrices products. After 
these kind of substitutions the output can be processed by MINCER. This procedure works 
well for all diagrams with no more than 2 gauge lines at operator vertex, i. e. we calculated 
668 from 682 three-loops diagrams^- Elapsed computing time was equal to 400 hours of 
CPU time. 

However, for the diagrams containing three gauge lines in operator vertex the described 
procedure requires too much computer time. Fortunately, taking into account diagram 
symmetry there are only 7 such graphs (some of them are presented in Fig. [T]). 

Each vertex of this type includes six different structures corresponding to all permuta- 
tions of gauge fields (see Appendix of Ref. |32j). Taking into account the color symmetry 
we have 15 different Lorentz structures for all 7 diagrams, containing such operator vertex. 
Each of these 15 structures can be obtained from the same combination of momenta and 
Lorentz indices. The difference between structures consists only in momenta labeling. A 
template was created for a chosen momenta combination and then other structures were 
obtained by means of proper substitution of momenta with the help of sed stream editor. 
The obtained data was used to declare a Table in FORM as the momenta are fixed for each 
structure. It turned out that the usage of Table allowed us to speed up by several orders 
of magnitude the process of generation of scalar expressions to be processed by MINCER. 
However, even after such an optimization the most complicated structures, which appear in 
the first three diagrams presented in Fig. [U required more than 30 hours of CPU time with 
FORM and more than 200 Gb HDD space (for temporary files). All diagrams containing the 
operator vertex with three gauge lines were calculated in approximately 250 hours. 

The final expression for 16th Mellin moment of the three-loop anomalous dimension of 

1 The results for all diagrams can be obtained from the authors upon request. 
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the flavor non-singlet transversity operator (pQ) was found to be given by 
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where a s = a s /(47r) = g 2 /(!6ir 2 ), £ 3 = £(3) is the Riemann zeta function, C A = N c and 
C F = (N 2 — 1)/(2N C ) are quadratic Casimir operators for gauge group SU(N C ), rij- is the 
number of active quark flavors and we put Tp = 1/2. This expression is in a full agreement 
with the prediction given in Ref . pL| , and can serve as a direct confirmation of the general 
result PP for arbitrary Mellin moment N. 

Acknowledgments 

This work is supported by RFBR grants 10-02-01338-a, 12-02-00412-a, RSGSS-65751.2010.2. 



5 



References 

[1] V. N. Velizhanin. larXiv: 1203. 10221 [hep-ph]. 

[2] A. K. Lenstra, H. W. Lenstra and L. Lovasz, Math. Ann. 261 (1982) 515. 

[3] J. P. Ralston and D. E. Soper, Nucl. Phys. B 152 (1979) 109. 

[4] R. L. Jaffe and X. D. Ji, Phys. Rev. Lett. 67 (1991) 552. 

[5] R. L. Jaffe and X. D. Ji, Nucl. Phys. B 375 (1992) 527. 

[6] J. L. Cortes, B. Pire and J. P. Ralston, Z. Phys. C 55 (1992) 409. 



[7] J. C. Collins, Nucl. Phys. B 396 (1993) 161 |arXiv:hep-p h/9208213|. 

[8] R. L. Jaffe and X. D. Ji, Phys. Rev. Lett. 71 (1993) 2547 |arXiv:hep-ph/9307329] . 



[9] R. D. Tangerman and P. J. Mulders, |arXiv:hep-p h / 9408305 , 
[10] D. Boer and P. J. Mulders, Phys. Rev. D 57 (1998) 5780 [arXiv:hep-ph/9711485| . 
[11] X. Artru and M. Mekhfi, Z. Phys. C 45 (1990) 669. 

[12] F. Baldracchini, N. S. Craigie, V. Roberto and M. Socolovsky, Fortsch. Phys. 30 
(1981) 505 [Fortsch. Phys. 29 (1981) 505]. 

[13] M. A. Shifman and M. I. Vysotsky, Nucl. Phys. B 186 (1981) 475. 

[14] A. P. Bukhvostov, G. V. Frolov, L. N. Lipatov and E. A. Kuraev, Nucl. Phys. B 258 
(1985) 601. 



[15] J. Blumlein, Eur. Phys. J. C 20 (2001) 683 [arXiv:hep-ph70 104099|. 

[16] A. Mukherjee and D. Chakrabarti, Phys. Lett. B 506 (2001) 283 
[arXiv:hep-ph/01020~03| . 

[17] A. Hayashigaki, Y. K anazawa and Y. Koike, Phys. Rev. D 56 (1997) 7350 
[arXiv:hep-ph/97072~08| . 

[18] S. Kumano and M. Miyama, Phys. Rev. D 56 (1997) 2504 |arXiv:hep-ph/9706420| . 

[19] W. Vogelsang, Phys. Rev. D 57 (1998) 1886 |arX iv:hep-ph/9706511] . 



[20] J. A. Gracey, Nucl. Phys. B 662 (2003) 247 |arXiv:hep-ph/0304113|. 



[21] J. A. Gracey, Nucl. Phys. B 667 (2003) 242 |arXiv:hep-ph/0306163|. 



[22] J. A. Gracey, JHEP 0610 (2006) 040 |arXiv:hep-ph/0609231|. 



6 



J. A. Gracey, Phys. Lett. B 643 (2006) 374 [arXiv:hep-ph/0611071] . 

J. Blumlein, S. Klein and B. Todtli, Phys. Rev. D 80 (2009) 094010 [ arXiv: 0909 .15471 
[hep-ph]]. 

J. Blumlein and J. A. M. Vermaseren, Phys. Lett. B 606 (2005) 130 
|arXiv:hep-ph/0411111] . 

S. A. Larin, F. V. Tkachov and J. A. M. Vermaseren, Phys. Lett. B 272 (1991) 121. 

S. A. Larin, T. van Ritbergen and J. A. M. Vermaseren, Nucl. Phys. B 427 (1994) 41. 

S. A. Larin and J. A. M. Vermaseren, Phys. Lett. B 303 (1993) 334 
|arXiv:hep-ph/9302208] . 

O. V. Tarasov and A. A. Vladimirov, Sov. J. Nucl. Phys. 25 (1977) 585 [Yad. Fiz. 25 
(1977) 1104]. 

A. A. Vladimirov, Theor. Math. Phys. 43 (1980) 417 [Teor. Mat. Fiz. 43 (1980) 210]. 

0. V. Tarasov and A. A. Vladimirov, JINR-E2-80-483. 

1. Bierenbaum, J. Blumlein and S. Klein, Nucl. Phys. B 820 (2009) 417 
jarXiv:0904.3563l [hep-ph]]. 

M. Tentyukov and J. Fleischer, Comput. Phys. Commun. 132 (2000) 124 
|arXiv:hep-ph/9904258] . 

P. Nogueira, J. Comput. Phys. 105 (1993) 279. 

T. van Ritbergen, A. N. Schellekens and J. A. M. Vermaseren, Int. J. Mod. Phys. A 
14 (1999) 41 larXiv:hep-ph/9802376] . 



J. A. M. Vermaseren, arXiv:math-ph/0010025, 



S. G. Gorishnii, S. A. Larin, L. R. Surguladze and F. V. Tkachov, Comput. Phys. 
Commun. 55 (1989) 381. 



7 



